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Abstract 

A model for the fragmentation of a nucleus is developed. Parallels of the de- 
scription of this process with other areas are shown which include Feynman's 
theory of the A transition in liquid Helium, Bose condensation, and Markov 
process models used in stochastic networks and polymer physics. These par- 
allels are used to generalize and further develop a previous exactly solvable 
model of nuclear fragmentation. An analysis of some experimental data is 
given. 
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I. INTRODUCTION 



The behavior of the distribution of fragments in an intermediate energy nuclear collision 
has attracted a great deal of interest. A previous set of papers [|l|-|3| proposed a statistical 
model for describing this behavior. Each possible fragmentation outcome is given a particular 
probability, and the resulting partition functions and ensemble averages are known to be 
exactly solvable. In this paper we extend that canonical ensemble model to allow more 
freedom in choosing the weight associated with a particular nuclear partition. Formulas for 
the ensemble averages and the partition function derived in earlier papers are generalized, 
and a recursive formula for the evaluation of the coefficients of the partition function is 
developed which is useful in computing low and high temperature behavior of the models. 

The balance of the paper is concerned with the application of these results to models 
of nuclear fragmentation, as well as a comparison of these results both to other models 
in physics, and to other models of fragmentation. An explicit parallel between this model 
and Feynman's approach [[| to the A transition in liquid Helium is proposed. The weight 
given to each possible cluster distribution in a canonical ensemble model is shown to be 
similar to that used by Feynman in the cycle class decomposition of the symmetric group. 
Moreover, the main parameter, called the tuning parameter x in ref. [Q, which contains 
the physical quantities associated with cluster formation, is shown to have a correspondence 
with a variable in Feynman's approach related to the cost function of moving a Helium atom 
from one location to another. The variable associated with this cost function is related to 
that part of the parameter x which has to do with internal excitations in a cluster, i.e., its 
internal partition function. Another model which has a similar structure as the canonical 
ensemble model is found in the physics of polymer formation. Polymer formation, when 
modeled by Markov processes, is a special case of the canonical ensemble models. Indeed, 
canonical ensemble models can be derived from a particular type of Markov process models, 
a fact which can be exploited when interpreting the underlying physics of canonical ensemble 
models. 
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Other models of fragmentation are easily compared to the canonical ensemble model. 
Models based on partitioning alone [0)0 are discussed briefly. They are also a result of 
assuming a certain weight is associated with each nuclear fragmentation outcome. However, 
the choice of weight is simpler than the models proposed here. Models based on percolation 
studies | fr6|jr8| also have some similar features. Markov process models ]5| as already noted 
are not only similar to canonical ensemble models, but rather are the exact same models, 
simply derived from a different, phenomeno logical, point of view. This is fortuitous, as an 
analogy with Markov process models can provide the basis for choosing particular canonical 
ensemble models for the study of nuclear fragmentation. A final model of fragmentation to 
compare the canonical ensemble model with is the generalized canonical ensemble models, 
of which canonical models are a special case. Canonical models have many advantages over 
this proposed generalization, but such generalized models may be useful in studying exotic 
fragmentation situations. 

A section of this paper is devoted to the discussion of the thermodynamic functions 
of canonical ensemble models. Since the canonical models are derived from a statistical 
mechanics assumption, it is appropriate to consider the computation of the typical thermo- 
dynamic functions. After deriving the appropriate formulas, they are applied to the case 
of an ideal Bose gas in d dimensions. The critical point, present for d > 2, is discovered 
by plotting the specific heat vs. the temperature. The connection between the zeros of 
the partition function and the critical behavior is also considered, as the recursion relations 
allow for the computation of the zeros for some non-trivial cases. The zeros for d = 2 and 
d = 4 are computed, and empirically appear to lie on simple arc-like curves. 

The extended canonical models are finally applied to the problem of nuclear fragmen- 
tation. After reviewing the general behavior of a number of models, which are found to 
vary widely in their fragmentation behavior, we focus on a small number of models which 
seem to be appropriate for an ensemble description of the fragmentation of ^gAu at 0.99 
GeV/amu. Several models give excellent results, but the statistics of the experimental data 
are not sufficient to distinguish a particular model from the others considered. One particu- 
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lar model, although a poor model of nuclear fragmentation, has very interesting properties 
and is analyzed further. In this model, the fragmentation distribution favors fragments of a 
particular size, with a Gaussian falloff in the distribution for larger and smaller fragments. 

The paper is organized as follows: Section [TI| develops the canonical model for fragmen- 
tation and discusses parallels of this model with the A transition in liquid Helium, Bose 



condensation and polymer physics. Section |TJ discusses a variety of other models for the 
fragmentation of a nucleus. Alternative partitioning models and percolation models are 
reviewed briefly. Markov process models are introduced, and are shown to give the same 



distribution as the canonical models of section |II A| when the detailed balance condition 
is satisfied. Additionally, generalized canonical models are introduced, which includes the 
canonical models of section |II A| as a special case. The thermodynamic properties of frag- 
mentation models are considered briefly in section [IV]. Various thermodynamic functions 
are derived, and as an application the computation of the specific heat of a finite Bose 
gas is given. The zeros of the partition function, which can indicate the presence of phase 
transitions, are also investigated for a few models. Section |V| discusses the behavior of the 
ensemble averages for various models. A comparison with some experimental data is made. 
Concluding remarks are in section [VT] and some additional calculations are included in an 
appendix. 



II. MODELS OF FRAGMENTATION AND PARTITIONING PHENOMENA 

In this section, we review an approach to fragmentation and partitioning phenomena 
based on the canonical ensemble of statistical mechanics. In this approach the fragmentation 
of a nucleus, or, in general, an object, is viewed in a statistical way with a weight given to 
each member in the ensemble of all possible distributions. Mean quantities, correlations, 
and fluctuations are obtained by averaging various expressions over the ensemble using this 
weight. The model considered is not limited to descriptions of nuclear fragmentation, and 
the rest of this section is devoted to introducing other areas in physics which have used a 
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similar type of description. Specifically, Feynman's description of the A transition in liquid 
Helium and an example in polymer physics are discussed. 



A. Exactly Solvable Canonical Models 

Exactly solvable canonical models, which can be used for the study of fragmentation and 
partitioning phenomena, were developed in a previous set of papers |]-|3j . Each partition or 
fragmentation is given the weight 



A\ -£ 1 fx k \ n 



^,xi~n- y (i) 

Qa{x) ^ n k \ \k ) 

where n = (n k )A — { n i, ■ ■ ■ , n A) is the partition vector for the fragmentation or partitioning 
of the A objects into n k clusters of size k, and x = (x k )A = ■ ■ ■ ,%a) is the parameter 
vector with x k characterizing the group or cluster of [size k. The partition vector must 
satisfy the constraint J2 k =i kn k = A and we denote the set of all partition vectors 



N a = jvn Y.kn k = A^ 



(2) 



The parameter vector contains the underlying physical quantities such as the temperature 
T and the volume V. The probability condition for PA{n, x) 

Y / P A (n,x) = l (3) 

HeN A 

determines the canonical ensemble partition function Za{x) = Qa{x)/A\ when eq. ([I]) is 
substituted into eq. (|3]). 

Previous papers dealt with two particular models. When all the x k = x, the x model of 
ref. JL], the partition function takes on a simple form 

Q A {x) = x(x + 1) • • • (x + A - 1) = T ^ x) A) (4) 

For the case Xi = xy, x k = x, k ^ 1, the xy model of ref. pfl, 

/ A\T(x + k) /r, 
x{y - 1)] (5) 



Qa(x) = E f 
fe=i V 



k Fix) 



Detailed studies |J show that the results of the x and xy models are quite similar for all 
cluster sizes k > 1. 

This paper considers in detail more general forms for Xk- Here, we explicitly show how to 

evaluate the partition function by simple recursive procedures. For convenience, we rewrite 

eq. ([I]) by making the following substitution 

x k (A,V,T) = x(A,V,T) 

k Pk U 

so that the dependence on the physical quantities is contained within a single parameter, x, 

and that the thermodynamic dependence and cluster size dependence are separable. Then 

the weight is given by 

A\ A 1 /-r\ nfe 

Qa{x;I3) k =i n k- \PkJ 

This is not an unreasonable constraint on the parameters, and is easily satisfied by many 

models. For example, a previous paper developed the result 

V \ a v k B T T ] 

x = — — — exp 



v (T) " { k B T s T + TqJ 

Pk = k (8) 

where T is the equilibrium temperature, V is the freeze out volume, and Vq{T) = h 3 / (27rm p fcsT 
is the quantum volume, with m v the mass of a nucleon. The a u is the coefficient in a sim- 
plified equation for the binding energy of a cluster of size k, E B = a v {k — 1). The e is the 
level density parameter related to the spacing of excited levels and To is a cutoff tempera- 
ture for internal excitations. In a Fermi gas model, Eo and the Fermi energy are related by 
e = Aep/ir 2 , and since Ep = p F /2m p can be obtained from 4(4:7ip F V/3h 3 ) = A, we find that 

k B T _ fn_\ 2 / 3 2m p k B T fV\ 2/3 

The evaluation of the partition function Za(x' : 0) = Qa(%] 0)/A\ and various ensemble 
averages from the weight can be derived from the generating function for Za{x; 0) 

oo f oo c I 

Z(n, = z A (%, P)^ A = exp £ J-n 11 (io) 



A=i 
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Using this function, it was shown in ref. || that the ensemble averaged cluster distribution 
(n k ) is given by 

/ \ x Z A _ k (x, (3) 
Pk Z A (x,P) 

More generally, if we define 

[z] k = z{z - 1) • • • (z - k + 1) (12) 

it was shown that 

fx \ kj ) Z A-J2. jkj (x-J) 

where Z k (x; (3) = for k < 0. 

The constraint J2k=i k{ n k) — A then leads to a simple recurrence relation for Za(x, (3) 

Z A (xJ) = ^J2Z A ^ k (xJ)^ (14) 

with Zq(x,(3) = 1. Then Zi(x,(3) = x/(3\, and so on. We can now calculate any ensemble 
average of n k using eqs. (0), (|i~4T) 

From the last equation we see that Z A (x, (3) is a polynomial in x of order A. To encourage 
this point of view, we will drop the dependence on (3 from the notation for Z A , making the 
dependence tacit. Then, the partition function can be written as 

Z A (x) = J2Z A k) x k (15) 
k=i 

where the coefficients Z A ^ can be determined from the recurrence relationship as follows. The 
first coefficient, Z^ A , is determined by the last term in the recurrence relation, (x/A) (AZq(x) / (3 A ) 
x/(3 A . So 

2? = 1 (16) 

Pa 

From this coefficient we can determine all the others 
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j=l Pj k=l 

tE E J-^V (i7) 

71 i=l k=2 Pj 



Thus all the coefficients Z A k ^ can also be obtained recursively 

zf = - A E (is) 
71 j=i Pj 

For fc near A, this expression can be used to obtain exact results for the coefficients. As- 
suming Pi = 1, which can always be done by redefining x, fi k such that x/ [3k is unchanged, 
i.e. x -> x//3 ± , f3 k -> /3 fc //?i , we find 

^) = J_ 
A! 

1 



7 l 



7 (A-1) 



7 (A-2) = - 

A " (A-3)!/3 3 2(A-4)!/3f 



(A-2)!ft 



1 1 



7 (A-3) 



1 1 



(A-4)!& (A-5)!&/3 3 6(A-6)!/3| 



7 (A-4) 1 , 1 / 1 , _J_\ /I OA 

A "" (A-5)!/% + (A-6)! \2Pi + m) [ } 

1 1 

+ ^ ^77^ + 
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2(A-7)!/3|/3 3 24(4-8)!$ 
1 1/1 1 



(A-6)\fo (A-7)\\p 5 {3 2 (3 4 (3 3i 
1(1 1 \ 



(A -8)! V2/3 4 /3| 2& 2 /? 2/ 
1 1 

+ 771 7 + 



6(A-9)!/M 120(A-10)!/3| 
In general we see Z^ 4 fc ^ depends on f3±, . . . , (3k+i- 

^"Sw-*-* E iKTW (20) 

n 



E>n ir = fc + s 



r=l 



The recurrence relation given by eq. (|14]) is simply solved for the case /3 k = k (as previ- 
ously noted) which gives 

Z A (x;(3 k = k) = ±J2\sW\x k (21) 

where S A axe Stirling numbers of the first kind. This model was analyzed extensively in 
refs. 0-0. 

Another case which reduces to a simple polynomial is /3fe = 1 which gives 

Z A (x;f3 k = l) = jxL 1 A _ 1 (-x) (22) 



with a Laguerre polynomial. The [3 k = 1 model is considered in detail in [12 



as a 



model for fragmentation and in |j] as an example of a Markov process model for clusterization 
of one dimensional objects. 

A final example whose coefficients are common mathematical functions is the choice 

/3 k = k\, 

Z A (x;P k = k\) = ^-J2S^x k (23) 
A - k=i 

with S A Stirling numbers of the second kind. This choice for the case x = 1 was considered 



in detail in ref. [0. It will also be analyzed more generally in section |V A . 

For any choice of /3 k , the recursion relation given in eq. (|14D holds. However, for some 
f3k there are simpler recursion relations. For example, if Q A (x) is given by an orthogonal 
polynomial, (e.g. /3 k — 1 ), then Q A+1 (x) = (a A + b A x)Q A (x) - c a Qa-i(x), as given in 
Abramowitz and Stegun 0. Table | lists some of these models. Note that the last choice 
for (3k in table [I] can be related to the Catalan numbers 

«-irr(?) 

Specifically, (3 k = 2 k ~ l / 'C k -i- If we use Stirling's approximation for the factorials in f3 k for 
this choice of /3 k , then (3 k ~ 2 k k 3 ^ 2 for large k. 

All the cases considered so far are special cases of some general forms The cases f3 k = 
1, k, k/2 k ~ l (^ {k zl ] ) 1 can be realized from 
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where [z] k is defined as 

[z] k = z{z + l)---{z + k-l) (26) 
The case fa = k\ is a special case of 

3 - fc![c]fc_1 (27) 

A - ^pr ( 27 ) 

One case which does not reduce to a commonly known polynomial and is of general 
interest is fa = k T . For example, an ideal Bose gas in d dimensions can be modeled by this 
choice, with r = 1 + d/2 (see section [IV B|) . In the large A limit, the Z A coefficients for 
small k are equal to z A / @Ai where the z A are only weakly dependent on A. Table [□] gives 

(2) 

the large A limit for these coefficients. Notice that lim^-soo^A = C( r )- 

Partition functions which satisfy recurrence relationships are common in physics. For 
example, consider the ideal Boltzmann gas 

Z A (x) = ^x A = jZa^x) (28) 

This is equivalent to the x model with the following choice of parameters 

V 

fa = S kl (29) 

Since fa = for k 1, this model has only fragments of size one, i.e., there is no clusteri- 
zation. For a noninteracting Boltzmann gas, this is clearly the correct behavior. 

Another example of a recurrence relation in statistical mechanics is the interacting Boltz- 
mann gas. For this example, Feynman [|J showed that when the three body and higher order 
terms are neglected, the spatial part of the recurrence relation 



Z A+1 = V\l--j Z A (30) 

where a = / °° — e~ v< " T ^ kBT ^j A7i:r 2 dr, with V(r) the two body potential. 

We now discuss some applications to other physical systems and illustrate the parallel 
with fragmentation phenomenon. 

10 



B. Parallel with Feynman's Approach to the A Transition in Liquid Helium 



In this section we give some of the results of Feynman's approach H] for the A transition 
which are relevant for the analogy to be discussed. Further details of the results quoted can 
be found in [11 and the references therein. 

The starting point is the partition function obtained by a path integral, given by eq. (11.52) 
in ref. |§ 

3N/2 



-F/k B T J_ ( 2nm'k B T 
N\\ h 2 



X 



Pes N 
x exp 



J d?Ri ■ ■ ■ d 3 R,Np(Ri, ■ ■ • , Rn) 

-^-E(^-^)) } (31) 

where N is the total number of Helium atoms, m' is the effective mass, Ri the coordinate of 
the i'th Helium atom p(R\, . . . , Rn) is the potential contribution and P is the permutation 
operator. A given permutation among the particles is illustrated in figure f[[ and can be 
visualized as the atoms being connected by a set of edges, the edges forming polygons 
(cycles) of various sizes (cycle lengths). After some algebra and approximations, the partition 
function is reduced to 

1 (Nh k \ nk 



e 



- F/kBT= E jjjr^ 



neN N V k J 

j m'd 2 k B T » \ 
X6XP i 2^~ ^ J ^ ' 

where hk = (c/k 3 ^ 2 + 1/N)l k , with / the number of nearest neighbors per lattice site for 
a particular choice of spatial discretization. This result is arrived at by a random walk 
argument. Specifically, starting at a given atom, there are l k random walks in fc-steps, and 
the fraction of these that are close (end up at the origin) is inversely proportional to the 
volume in which the random walk is likely to end. This volume is in turn proportional to 
k 3 / 2 . The random walk result is then corrected for the very large polygons that encompass 
a large fraction of the sites. This determines the form for hk, as given above. 
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Comparing the results of this section with those of section |II A|, we see a strong parallel 



between the model of fragmentation and the model for the A transition. First, the exponent 
of eq. fl32|) is identical to that internal excitation function of eqs. (||), (|9]), up to a numerical 
constant, if we neglect the cutoff temperature factor (T — > oo) and make substitutions 
for the density (A/V — > d~ 3 ) and the quantum volume (v (T) — > d 3 ). Secondly hj, 1 is 
analogous to the parameter f3k in eq. (0). In fact, the modification Feynman makes to hk 
can be well motivated in the case of nuclear fragmentation, and we will consider the case 



1/Ac = a/k T + (1 — a)/k in section [V B| . Third, the partition function has a formal structure 
identical to the xy model. 

From the above remarks we note two important issues in the choice of weight given to 
each partition, fragmentation or grouping. One issue is the choice of j3k and the second is 
the relation of x to the physical quantities. 

C. Polymerization Processes 

In organic chemistry, one can model the formation of polymers using the same models 
considered here, if one is not overly concerned with the details of the molecules formed, only 
in their size. We assume that the molecules combine to form polymers by forming bonds 
between molecules, up to a maximum of / bonds on a single molecule. 

If existing bonds between molecules break at a rate k and new bonds form at rate 
proportional to the number of sites available for new bonds, then one can show (see ref . || , 
chapter 8, and references therein) that the equilibrium polymer distribution is given by 
applying the x model with the following parameters 

X = K 

_ jfc!((/-2)fc + 2)! 
Pk ~ ((/-l)*)! (33) 

In this case A is the number of molecules, (rife) is the expected number of polymers containing 
k molecules. Kelly obtained this distribution for polymer sizes by developing a Markov 
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process model for polymerization. We will show later in this paper (Section [III C|) that such 
models are equivalent to the x model with an appropriate choice for (3k- 
The Stirling limit of (3k for / = 3 is 

(3 k « ^2- 2k k 1/2 {k + l)(Jfe + 2) « ^2- 2k k^ 2 (34) 

Thus this model of polymerization is asymptotically similar to the Bose gas choice for (3k for 
d = 3 for large k (the term 2~ 2k does not affect the distribution and can be ignored). Since 
the replacement (k + l)(k + 2) ~ k 2 is not very good for small k, this approximation is poor 
for large x, since then the configurations are dominated by small clusters, whose behavior 
is determined by (3k for small k. Taking (3k = k x l 2 {k + l)(k + 2) and x = k/\/tt, figure |2| 
shows a comparison of the mean number of polymers with k units with the original (3k and 
x of eq. (H) for / = 3,x = 40. 

III. COMPARISON WITH OTHER MODELS OF FRAGMENTATION 

In this section, we discuss a number of other models of nuclear fragmentation. Models 



based on partitioning alone ||i4| , |15f are similar to the model outlined in section [II A| but with 



a simpler choice for the partition weight. Percolation models [|T6| , P~8|1 are derived from far 
different assumptions. Markov process models are identical to the x model, but derived from 
a phenomenological point of view. They are useful for considering what forms of (3k would be 
appropriate for modeling fragmentation phenomenon. Lastly, a generalized canonical model 
is introduced. Although it does not allow for easy computation of the various ensemble 
averages, it may be useful in investigating the behavior of exotic partition functions. 

A. Models Based on Partitioning Alone 



Sobotka and Moretto |TJj] and Aichelin and Hufner (jOJ discussed a model of fragmen- 
tation based on partitioning alone with no microstate counting factor or tuning parameter. 
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In particular, their assumption is that every partition is equally likely and that alone de- 
termines the fragmentation. The number of partitions of A is P(A) which can be obtained 
from the generating function 

OO OO 1 

E P(A)x A = n -—a ( 35 ) 

A=0 A=l 1 x 

This is asymptotically given by the Hardy- Ramanjuan result 

P(A) « — L_ e^v^ (36) 

The number of partitions of A with fixed multiplicity m = J2k=i n k is m ) an d is given 
by the recurrence relation 

P(A,m) = P(A- l,m- 1) +P(A-m,m) (37) 

In this approach the frequency of clusters of size k is (nk) = P(A — k) + P(A — 2k) + . . . 
which can be reduced to 

(n k ) » (38) 

* ; exp{(^)V2fc} _l 1 J 

The above simple model of fragmentation can be generalized to include a tuning parame- 
ter x which contains the underlying physical quantities such as volume, temperature, binding 
and excitation energy associated with a fragmentation process. In this generalization of the 



models given in ref. fT^JTq] , the weight given to any partition is 

x m 

QAx) (39) 

Here the Qa(x) is the normalization factor for this model, and is given by 

Q A (x)= J2P(A,m)x m (40) 

m=l 

Once Qa{x) is obtained, various mean quantities can be found. For example 

1 A 

^ = TTua £ x r Q A -rk{x) (41) 
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and for k ^ j 



while for k = j 



(rijUk) = x r+s Q A -rk-sj(x) (42) 



(n k {n k - 1)) = -, : x r+s Q A -{r+s)k{.x) (43) 



For large A and xA ^> 1, (n fc ) can be shown to approach 

K) ~^p{(€f) 1/2 ^}-i (44) 

At x = 1, the above formula reduces to the result of eq. fl3~8|), as expected. 



B. Percolation Models 

The x model has one variable that describes the degree of fragmentation. The x ranges 
from to oo as the temperature changes over the same range. Another approach to cluster 
distributions is based on percolation which also uses one parameter. The application of per- 
colation to nuclear fragmentation was developed by several groups [I6,|I5|] . The percolation 



models are of two types: bond and site. The bond type assign a certain probability p of 
having a bond between lattice sites, while the site type assign a certain probability of having 
a site occupied. The x of eq. (H) and x/A have terms which deal with volume or density 
effects, and binding effects. Also included in x are thermal effects through vq and internal 
excitation energy considerations not present in percolation studies. 
In percolation studies, the number of clusters of size k is given by 

(n k (p)} = ^f((p- Pc )kn (45) 

where p c is the critical probability above which an infinite cluster exists, / is a scaling 
function, and r, o are critical exponents. 

For the choice (3k = k, the x model has a cluster distribution approximately given by 
(see ref. @) 
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(n fe (A;«Ae«l))^7 e " efeM ( 46 ) 

k 

with x = 1 + e, giving t = 1, a = 1. At e = or x — 1, n*, = 1/k, a hyperbolic power law 
behavior. Other choices for would give a different k dependence for n k . 



C. Markov Process Models 

An alternative point of view for modeling nuclear fragmentation comes from Markov 
processes, which allows the underlying physical phenomena to be reflected in the equilibrium 
distribution. The idea is to consider a method by which a cluster configuration can change 
into another configuration, and then to derive what the equilibrium distribution is for such 
a method applied to the set of states. Rather than assuming what the probability for 
each state is, it is derived from the distribution achieved by applying the Markov process 
repeatedly. 

For example, we can consider that the underlying physical processes are the joining of 
two fragments to form a new larger fragment and the splitting of a larger fragment into two 
smaller fragments. Of course, fragments joining and breaking into more than two groups 
are possible, but we will ignore that for now, assuming that those processes are of lesser 
importance and will not materially affect the overall equilibrium distribution. We denote 
these processes by a transition operator T which acts as follows on n, the configuration 
vector. 



Ttn = 




n 2 , ■ . 


.,Uj 


-1,.. 


. ,n k - 1, . 


• • , ni + 


1,. 


■■) 


Tpn = 




n 2 , . . 


.,Uj 


-2,.. 


• ,ni + 1, . 


••) 






T] k n = 


(ni; 


n 2 , . . 


.,Uj 


+ 1,.. 


. ,n k + 1, . 


• • , ni - 


1,- 


••) 


T] 3 n = 


(ni; 


n 2 , . . 




+ 2,.. 


.,ni - 1,. 


••) 







Now obviously these operators do not conserve particle number unless j '• + k = I, so we need 
only consider a smaller set of operators 
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j-\~k 

T 3k = T$ k (48) 

Suppose these processes occur at some rate, denoted by q(n, n') for n transforms to n' . 
For example, 

q(n,T jk n) = X jk nj(n k - 5 jk ) 

q(n,T jk n) = fx jk n j+k (49) 

This is a very reasonable choice, as the probability is proportional to the number of fragments 
available for such moves. If there are none, then the transition probability is zero, as needed. 

We expect that this process when applied repeatedly to a configuration will lead to 
an equilibrium configuration, at which point the rate at which transitions occur to a new 
state, weighed to reflect the equilibrium distribution of the original state, is equal to the 
rate at which transitions occur back to the original state, weighed to reflect the equilibrium 
distribution of the new state. In other words, if Px(n) is the equilibrium distribution, it 
must satisfy the detailed balance condition: 

P A (n)q(n, T jk n) = P A (T jk n)q(T jk n, n) (50) 

If there exist positive numbers ci, . . . , c A such that 

CjC k Xj k = Cj +k fij k (51) 

then it is easy to show that the equilibrium distribution is given by 

1 n k 

PA(n) = ^U C -h (52) 
^ a k n k \ 



If for some choice of Xj k , fij k , we get c k = x//3 k , then the models considered in section |II A 
are reproduced. Once we know c k we can use eqs. ([13]), (0) to solve for the various ensemble 
averages. In fact, the recursion relationship stated above is another way of expressing the 



recursion relationship relating partition functions developed in section [11 A| . Of course we 
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still have not produced a set of Xjk, fijk which satisfy eq. ([5ID . A very general solution 
(though not unique) is given by the following choice 

X 3 k = af(j)f(k) 

Vjk = Pf(J + k) (53) 

where f(k) is any nonnegative function of k. It can be easily shown that the solution to this 
model (in the language of the x model) is 

x = (3/a 

Pk = /(*) (54) 

Thus we can always generate the x model from a Markov process, using the above 
prescription for Xjk, fijk- The benefit of this approach is that the connection between the 
underlying physical processes and the corresponding probability distribution is easier to 
consider. The parameters Xjk determine the association aspects of the model, while the fijk 
determine the dissociation aspects. Table fTJ lists a number of examples and their solution 
in terms of the x model. 

Suppose we think of the fragmentation process as having all the nucleons involved con- 
strained to move in a small volume of space for a period of time long enough for the fragments 
to achieve an equilibrium. The probability of joining two fragments should be determined 
mostly by the density of the nucleons in this volume, and the cross section of each fragment. 
We expect larger nuclei to accrete smaller nuclei due simply to their larger cross section. 
Therefore Xjk should increase monotonically in j and k. All the nuclei created can break 
up into smaller nuclei. We expect larger nuclei to be more unstable than smaller nuclei. 
This is not strictly correct, as larger nuclei are energetically favorable in their ground state 
configurations. However, in the aftermath of a high energy collision, the added excitation 
energy and angular momentum should make larger structures unstable. Therefore Hjk should 
increase monotonically in j and k as well. This suggests that the above example might be 
a reasonable model of nuclear fragmentation, provided f(k) increases monotonically. 

18 



Since all the transition operators can be generated from a smaller set of operators, an 
obvious simplification of this model is to consider the distributions generated by a restricted 
set of transition operators. For example, we can restrict the Markov transition operators 
down to two operators, one which adds a single nucleon to a fragment, the other which 
removes a single nucleon from a fragment. We denote these restricted operators by 

T k = 7\, fe _ a 

T k ee r 1 '* -1 (55) 
and the restricted transition coefficients 

fik = j"i,fc-i (56) 

Using these restricted operators any general operator (and therefore any state) can be gener- 
ated by composing the restricted operators together in an appropriate combination. Specif- 
ically, 

T jk = II T j+k+2 ^ f[T m f[T n 

1=2 m=2 n=2 

j k j+k 

T jk = 11 T j+2 ~ m JJ T k+2 ~ n JJ Ti (57) 

m=2 n=2 1=2 

Therefore if the recursion relationship given by eq. ([51]) holds, it can be built up from 

Cfc/ifc = c fc _iCiA fc _i (58) 

For these models, one needs only specify the parameters X k and The full Xjk, l^jk can 
then generated from the set of c k by applying the full recurrence relation eq. (|5lD . A number 
of models of this type are listed in table 0. 

Although these models give the same ensemble averages as the more general models, they 
have different underlying physics. The parameters fi k , A& still determine the dissociation and 
association aspects of the model; however it is easier to discover functions that satisfy the 
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restricted recurrence relation given by eq. (|58|) than it is to satisfy the full recurrence relation 
of eq. fl5"T|). In other words, these models allow more freedom to adjust the dissociation and 
association rates. For instance the model 



Hh = Pk T (59) 



which has the solution 



x = — 
a 

is not easily produced from an obvious choice of parameters Xjj., 

The above model conforms to the idea introduced earlier of the nucleons combining 
and dissociating as a gas in a limited volume. Single nucleons can detach themselves from 
larger nuclei, and larger nuclei can accrete single nucleons. If we assume that the nucleons 
automatically combine if they get close enough, then the association rates should be roughly 
proportional to the surface area of the large fragment, i.e a = 2/3. If we assume that nucleons 
near the surface of the large fragments leave at a constant rate, then the dissociation rates 
should also be a surface phenomenon, and r = 2/3, which implies (3k = k 2 ^ 3 . We will see 
in section |VB| that this model does not reproduce experimental data well. Instead models 



with a = r > 1 produce better results. 



D. Generalized Canonical Models 

As a final example of a model of fragmentation, we discard the notion that the model 
must be derived from any particular choice of weight. Indeed, the usual process of choosing 
a weight and then deriving its partition function can be reversed. A partition function can 
be chosen, and a weight scheme that generates such a partition function can be computed. 
It is important to note that the choice of a partition function is not sufficient for fixing such 
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a weight scheme. Many different choices for a weight lead to the same partition function 
and only additional assumptions can fix the weight scheme. For example, the model given 
by the weight 



PA(n,x) = ——ll — (- 
Qa{x) £*i n k \\k 



and by the weight 



P A (n,x\m = Y^ n-k) 



x' 



Qa(x) Sf 



both give the same partition function 



(61) 



(62) 



Qa(x) = £ fi£ 



(m) 



m=l 



X 



Suppose that we are given a partition function 



QA{x) = Y.Q7x k = Y,W A {n) 

k=l {n} 



(63) 



(64) 



and we want to determine a weight scheme that generates this partition function. For the 
canonical model with (3^ = k, the weight is given by WA_{n) = M 2 (n)x m , where m = J2 n k, 
M 2 (n) = A\j Y[nk\k nk . An obvious generalization of this weight would be the choice 

A 

k \n k \ \k) 

i.e., WA(n) is equal to the standard canonical model weight, up to a factor that depends 
only on the number of fragments. For this weight we can easily show 



W A (n G N A \m = ^> fe ) = W A m) A\ f[ 



(65) 



E W A {n) = W$ 



(m) 



? (m) 



X 



Q { fx m 



(66) 



which implies that 



or that the weight is given by 



W 



(m) 



s 



(m) 



(67) 
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W A (n E N A | m = ^ n k 



rt(m) 



71 fc 



(68) 



Now if we could write down the generating function for Qa{%) we could compute (rife) 
as was done in section [II A| . However, in general there is no generating function for Qa(x). 
We can calculate (m) = J2( n k) exactly, however. 

x dQ A 



(m) 



(69) 



Qa(x) dx 

The (rik) results are not entirely inaccessible. They can be obtained by a Monte Carlo 
simulation of the partition function. Note first that the partition function can be expressed 
as a sum over the permutation group. 



Qa(x) 



E 



A\ Q ( a ] ^ 



r, ev ,v X m=in k \k^ \sl 



M 



x 



E 



X 



m(p) 



E e 



-s(p) 



(70) 



with 



5(p) = log 



5 



(m(p)) 



logQ^ Tl ^' ) — m(p) logx 



(71) 



We can simulate this action over the set of permutation using the Metropolis algorithm, 
with any of a number of choices for transition functions. For example, if we denote p k as 
the action of the permutation operator on k, one possible transition function is p k <-> Pk+i 
for a uniformly random choice of k. 

The choice of Qa{%) cannot be completely arbitrary. In the large x limit, it must reduce 
to Qa{x) = x A + 0(x j4 ~ 1 ), which fixes = 1- Similarly, the small x limit fixes the choice 
of Qa ■ Other than these considerations and the requirement that that the coefficients be 
positive, there are no other restrictions on the partition function. Indeed, some fairly exotic 
choices can be made. 

One example of a generalized canonical model is given by the partition functions gener- 
ated by 
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Qa+i(x) = Q A (^-^)(cx + d) A+1 Q A (72) 

with Qi(x) = x, and 9^ a constant which corrects • Models of this type satisfy a 
simple recurrence relation, but that recurrence relation is quite different than the canonical 



recurrence relation given in eq. (H). Another model, even more exotic, is given by 

Q2a{x) = Q A { qX{X ^ l) ){x + <1? A ®a (73) 
x + q 



Models such as this are interesting in studies of the roots of partition functions p5| , p6| , since 
the computation of large numbers of zeros for such partition functions is easily accomplished. 
The distribution of roots in the complex plane can have a fractal character when the re- 
quirement that the coefficients be positive is relaxed. The case q — — 1, shown in figure |3|, 
reveals one such complex distribution of roots. The roots lie on the boundary of a series of 
copies of the Mandelbrot set. A discussion of zeros of the partition function for canonical 



models will be made in section IV C 



IV. THERMODYNAMIC PROPERTIES OF FRAGMENTATION MODELS 

In this section we consider the various thermodynamic functions that can be computed 
from the partition functions obtained from the weight given in eq. (|7|). As an example, the 
specific heat for a finite Bose gas is computed. Finally, a discussion of phase transitions 
leads us to consider the location of the roots of the partition function on the complex plane 
for several models. 

A. Thermodynamic Functions 



In section |II A| we introduced a model with the thermodynamic variable confined to a 
single parameter x. Making this assumption allows us now to simply calculate the thermo- 
dynamic functions of such partition functions. Since 

e -F { A,V,T )/kB T = Qa{x) = Al ^ JJ J_ * (74) 

{nk}k=1 n k l \Pfc/ 
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is the partition function for a thermodynamic system with x = x(A, V,T), it is straightfor- 
ward to calculate the thermodynamic functions from the free energy. First, let us introduce 
dimensionless variables for T and V 

T V 

T\ and V\ are arbitrary reference points, but we will find it convenient (in nuclear fragmen- 
tation) to use the values 



where a v defines the scale of binding energies, and r is the classical radius of a nucleon. 

We can express the various thermodynamic functions in terms of (to), (to 2 ), x and its 
derivatives. The calculations are simple, and here we quote the results for the dimensionless 
entropy, pressure, free energy, energy, and specific heat. 



k B Ti = 



a v , Vi = 



^r 3 A 



(76) 




(77) 



(78) 



t\nQ A (x) 



(79) 



1 



(F + TS) = f + ts 



u 



l — 



k B T x 




(81) 



(80) 
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where the ensemble averages of m, m are obtained as follows 



<H?*) = «OT (82) 



mm 



i)> = (e»*(x>*-i) 



X2 ^ A (83) 



Qa(^) dx 2 



B. Specific Heat of a Bose Gas 

We can apply the above expressions to obtain the specific heat of a finite Bose gas. A 
Bose gas in d dimensions can be modeled by the x model, with 

V 

x 



vo(T) 

Pk = k 1+d ' 2 (84) 

where v (T) = (ti 2 /2Trni p k B T) d / 2 . Combining eqs. (j8l|) , (|8lD we arrive at a simple formula 
for the specific heat 

c v = ±(m} + ^((m 2 }-(m} 2 ) (85) 

For d > 2, there is a phase transition (Bose-Einstein condensation) in the infinite particle 
limit, which can be seen as a cusp in the specific heat at the critical point, x c = A/((d/2). 
For finite gases, the partition function is smooth, so there is no cusp. However the specific 
heat does reach a maximum near the critical point, suggesting the cusp will onset in the the 
large particle limit. This behavior is illustrated in figure f|. 

This model could be taken as a model of nuclear fragmentation, with a different expres- 



sion for x. Indeed, it will be shown in section [VB| that models with (3k = k T with r > 1 
are fairly good models of fragmentation. The analog of Bose-Einstein condensation into 
the ground state is the formation of the largest cluster, that is the "condensation" of the 
nucleons into the single cluster with k = A, which is called the fused mode in ref. fL|. 
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The canonical ensemble partition function is obtained from eq. (|14]) . The polynomial 
associated with this partition function, given by eq. (|3~5|) , has the following coefficients which 
are obtained by applying eqs. flT6|), (|T8|), and (|20|). 



(,4) 



^4 



1 

~A\ 



2 i+rf/2(^_ 2 )! 



7 (^-2) _ 1 | I (o 6 \ 

A - 3i+d/2( A _3), -r 2(2 1 + rf / 2 ) 2 (A-4)! 1 J 



7 (A-3) 



1 1 



4i+d/2(^_4)i 2 1 + af / 2 3 1+d / 2 (A - 5)! 
1 



Z (2) 



6(2 1 + d / 2 ) 3 (A-6)! 

C(l + d/2) 



^ A 



1 



^l+d/2 

These coefficients can be used to obtain the behavior of various thermodynamic functions 
in the low and high temperature limits. 

In the following subsection, we will consider another way of obtaining the phase transi- 
tions of a partition function. Using the results obtained in earlier sections, we attempt to 
calculate the zeros of the partition function for various A. 



C. Zeros of the Partition Function and Phase Transitions 

The canonical partition function Za{x) is a polynomial of order A in x with positive 
coefficients. For example, for /3k = k, Z&(x) = x(x + 1) ■ ■ ■ (x + A — 1)/A\ which when 
expanded gives 



1 A 



m=l 

,(m) 



Z a{x) = ^ E l^r * m (87) 



with S A Stirling numbers of the first kind. By a theorem of Gauss, a polynomial of order 
A has A roots or zeroes in the complex plane. For positive coefficients, no real roots are on 
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the positive real-axis, which is also the physical meaningful axis. The above example has its 
roots at the x = and the negative integers x = — 1, —2, . . . , —A + 1. 

Complex roots correspond to an extension of the real temperature into the complex 
plane. Lee and Yang, in their discussion of phase transitions, showed that such transitions 
manifest themselves as zeros of the partition function approach the real positive axis as the 
thermodynamic limit A — > oo is approached. Taking the logarithm of the partition function 
to obtain the free energy can then lead to a singularity. 

To illustrate the above remarks, we consider the example of the 2-d Ising model on an 



m x n lattice. Kauffman (pjj showed that the partition function for this model is given by 



where 

f rs = cos(27rr / m) + cos(27rs/ra) 
v = tanh(J/k B T) (89) 

In this case the zeroes of the partition function are located on two circles in the complex v- 
plane, namely v = ±1 + \/2e lb '. The physically meaningful domain of the v -plane is the part 
of the real line v G (0, 1) ( assuming J > 0). The zeros of the partition function approach 
this domain at one point, v = (—1 + y/2). This implies there should be a phase transition 
when 

Wc = r^- =- = 2J i — r (90) 

tanh- x (-l + y/2) log(l + v^) 

which is the commonly known value for the critical temperature. 

Now consider the zeros of the x model partition function in the complex x-plane. The 
physically meaningful domain of x, the part with positive temperature, is the positive real 
axis. So for a phase transition to manifest itself in the infinite particle limit, the zeros of 
the partition function must approach the positive real axis. The fa = k model therefore has 
no phase transition, for in the infinite particle limit, the roots of the partition function are 
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zero and the negative integers, which never approach the real temperature domain. Another 
example we consider is (3k = k T for the cases r = 2, 3. This corresponds to an ideal Bose 
gas in two and four dimensions by eq. fl34"|). So we expect that the zeros should approach 
the physically meaningful domain for large A for the case r = 3, but not for the case r = 2. 
Figures |5| and ^ illustrate the zeros for the models (3k = k 2 and (3k = k 3 for A = 25, 50, 75, 100. 
These graphs suggest that for both models, the roots lie on simple curves. Whether these 
curves close on the positive real axis is not clear from these small A results. The roots near 
the negative real axis scale with A, which suggests that the crossing point on the positive 
real axis should scale with A as well. This agrees with the known behavior of the critical 
point, which is given by x c /A = l/((d/2). The fact that the curve will not close for d = 2 
is not evident from the graph, however. 

V. GENERAL BEHAVIOR OF MODELS AND COMPARISON WITH 

EXPERIMENTAL DATA 

In this section we consider the general behavior of (rife) for various (3k and how well 
these models proposed actually fit some experimental data obtained from heavy ion collider 
experiments ||28|| . This data was obtained from emulsion experiments for ^gAu at 0.99 



GeV/amu. 415 events were recorded and identified by the charges of the fragments, i.e. 
each event is represented by n = (m, . . . , njg) where n z represents the number of fragments 
with charge z in a given event. Ensemble averages were obtained by averaging over all 
events. 

Since the experimental method only measured the electric charge of fragments leaving the 
collision, there is some question about how applicable are models developed considering only 
nucleons with no separation into protons and neutrons. Since the nuclear force treats protons 
and neutrons nearly identically, and the models proposed derived from a combinatorial 
viewpoint, the models should be identical whether one includes the neutrons or not. In fact, 
it can be shown that for the simple f3k = k model, that the results are the same whether 
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one considers Z nucleons coming out, or whether one considers A nucleons coming out, 
but only the Z protons can be followed, such that one must sum over the possible neutron 
configurations to obtain the expectation values. 



A. General Behavior of the Models 

Because all the models must satisfy J2k kn k — A there are restrictions on the form of 
the distribution. If we graph (mk) = k(nk) vs. k, then the area under the curve must be 
equal to A. For different choices of x and fik, the area will be distributed differently. In this 
subsection we discuss typical distributions for these models. 

All x models have simple behaviors which are easy to obtain at large and small x. For 
small x, all models will produce (m^) with most of the area under (m^). This is because 
for small x, the partition function is given almost entirely by the term proportional to x. 
This implies that one fragment is the most likely outcome. For k ^ A, (n^) and (m^) are 
proportional to x since 

W=^«^ (91) 

Pk ^ A PkPA-k 

For large x, all models will produce (m^) with most of the area under (mi). This is because 
for large x, the partition function is given almost entirely by the term proportional to x A . 
So A fragments is the most likely outcome. For k ^ 1, (n^) and (mj.) are proportional to 



x 1 k since 



_ x Z A _ k (x) x 1 k A\ 
{nk) ~(3 k Z A * (3 k (A-k)\ {J2} 

The model f3k = k was discussed in an earlier set of papers. For small x, most of the 
area is near k = A, as expected. For x < 1, (m^) is monotonically increasing. At x — 1, 
(mk) = 1 for all k. For x > 1, (rrik) is monotonically decreasing. At large x, most of the 
area is near k — 1, as expected. This is shown in figure 0. 

For models with ftk = I, for i C 1 most of the area is below k = A. As x increases, 
some area is distributed along the rest of the graph, mostly around k = A/2. As x keeps 
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increasing, the area continues to be redistributed, until most of it is distributed about a 
point k < A/2. At large x it takes on the usual distribution. This is illustrated in figure 

For the model (3k = k T , with < r < 1, the behavior is very similar to = 1. For 
small x, (ink) increases monotonically,with most of the area near k = A. As x increases, the 
right hand side is diminished till the graph attains two turning points, a local minimum near 
k = A, and a local maximum at a point k < A/2. The local minimum soon disappears, and 
the local maximum migrates left till it reaches k = 1 at large x. This behavior is illustrated 
in figure ||. 

The models (3k = k T , r > 1 are all very similar. For small x, (rrik) starts monotonically 
decreasing, reaches a minimum, then near k = A rises rapidly. As x is increased, the 
small k behavior is unchanged, but eventually the large k part turns downward. Thus for a 
small range of x, the graphs have two turning points. As x continues to increase, the local 
maximum eventually disappears and the typical large x behavior onsets. This behavior is 
shown in figure |T0| . 

The model (3k = k\, for (wifc) is mostly under k = A, as expected. From eq. (PT|). 

we see (n^) ~ x (a!) f° r k A, so the remainder of the area is binomially distributed 
about k = A/2. As x increases, the amount of area under k = A diminishes, the balance 
appearing around k m A/ 2 in a binomial or Gaussian distribution. Once most of the area 
has disappeared from k = A, the Gaussian distribution at the center moves to the left as x 
is increased. For very large x, most of the area is under k = 1, as expected. 

A simplified description of the (3k = k\ model can be obtained by making the following 
approximations. For a given x, the partition function is strongly peaked about a particular 
number of fragments m (i.e. Qa{x) ~ Q^x m ). To estimate this value of m we first use the 
approximation for the Stirling numbers of the 2nd kind, S™ ~ m A /m\, in the expression 
for the partition function. This approximation is very good for m < A, and is only off by 
a factor of two for m ~ A/2. Using this approximation, we can calculate which term in the 
partition function dominates, i.e. for what value of m , Q^ 1 x m ° is maximal. We discover 
the following nonlinear equation for itlq 
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m = xe A/mo = m {A, x) (93) 

In obtaining this result, Stirling's approximation for m\ ~ m m+l ' 2 e~ m y2m was used. The 
distribution (rik) is obtained from eq. (|ll|)- If we niake the approximation 

x Z A ^ k (x) (A\ _ k 

{nt) = e.^W * [k) m * Mx) (94) 

where = mo (A — k,x). We now make the assumption that 

rrik ~ mo(A — (k),x) = m (95) 

The value of /a(x) can be obtained by imposing k(nk) = A, which gives 

(n k ) w ^p fc (l-p) A - fe (l + m) (96) 

where p — l/(m+ 1). This distribution is binomial, with 

(A;) = (97) 
w m + 1 

Using this result we can determine the equation for rh. 

m = m (A- , x) = m (A^—, x) = xe A/ ^ +l) (98) 

m + 1 m + 1 

Figure [TT] compares the exact behavior with this approximation. From the figure we see 
that (nk) is reasonably well described by this approximation. 



B. Experimental Comparisons 

Nuclear fragmentation has a characteristic distribution which is met generically by only 
a few of the above models. For the experiment we will analyze, {rrik) drops, then rises. This 
suggests models with fa = k T with r > 1 might be satisfactory if one x is used. Models with 
two or more x's are considered in ref. 0. The large k behavior is somewhat indeterminate. It 
could be rising or falling; there are not enough events to determine the behavior accurately. 

Fits were made to log(n&), dropping from the experimental distribution any (rik) that 
were zero due to insufficient statistics. A previous paper || showed that for x = 0.3, the 
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(3k = k model gives a fairly good fit. A better fit is obtained by using two or more x's. 
Here, we consider the models proposed in the previous subsection, as well as a mixed model 
analogous to Feynman's choice for the A transition [|J, l//3k = a/k+(l — a)/k T , with various 
t. These results are shown in table M, and in figure [12 



VI. CONCLUSIONS AND SUMMARY 

This paper presents a detailed investigation of a set of exactly solvable canonical ensem- 
ble models of fragmentation processes and discusses some of its parallels with other areas. 
Specifically, parallels between the description of the fragmentation process and other areas 
are developed which include Feynman's approach to the A transition in liquid helium, Bose 
condensation, and Markov process models used in stochastic networks and polymer physics. 

The partition functions derived from various weights given to each member of the canon- 
ical ensemble, are shown to be polynomials in a parameter x. Simple recurrence procedures 
are developed for obtaining the partition function and the coefficients in the associated 
polynomials. The variable x, called a tuning parameter, contains the underlying physical 
quantities associated with the description of the the different processes considered. For ex- 
ample, for the ideal Boltzmann gas, ideal Fermi-Dirac gas, and ideal Bose gas, x involves 
the thermodynamic variables V (volume) and T (temperature) through the quantum volume 
vo(T). For fragmentation, x also includes binding energy and excitation energy coefficients 
associated with cluster formation and in the Feynman description of the A transition, x 
contains the cost function of moving a helium atom from one position to another. This cost 
function for the A transition is shown to be related to that part of x in fragmentation pro- 
cesses that involves internal excitations in a cluster. The length of the cycle of a permutation 
in the symmetrization of the Bose system wavefunction is the analog of the cluster size k 
and the cycle class decomposition of the symmetric group is equivalent to the partitioning 
or grouping of the original A objects into clusters of various sizes. 

Besides the tuning parameter x, the weight given to each member of the ensemble con- 
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tains a quantity called (3k which gives the cluster size or cycle length k dependence of this 
weight. Various choices for (3k are considered, and a wide range of different types of behavior 
can be found for different choices for (3k- The Bose gas in ^-dimensions has (3k = k l+d l 2 and 
Bose condensation exists for d > 2. The Feynman approach to the A transition is based 
on (3(( x = ak~ h l 2 + (1 — a)k~ l . A previous model of fragmentation [|T]-[3] used (3k = k, a 
choice leading to very simple results due to some theorems in combinatorial analysis. The 
partition function for this last choice is a simple polynomial in x whose coefficients are the 
signless Stirling numbers of the 1st kind. A more general form, (3k = k T , is also investigated 
here. For a partition function that leads to a phase transition, r > 2. Also considered here 
is the choice (3k = k\, which is shown to have some interesting properties. The partition 
function 'in this case involves Stirling numbers of the second kind, and the distribution of 
fragments obtained from this partition function represents systems which split into large 
equal or nearly equal size pieces. Specifically, this choice gives rise to Brownian or binomial 
type distributions of clusters whose peak is centered around cluster sizes A/m where A is 
the number of objects and m is the mean multiplicity, which is a function of x. 

A consideration of the thermodynamics of fragmentation systems led to an investigation 
of the behavior of the partition function when x is a complex number. In particular, the 
zeros of the partition function are studied in the complex plane and the connection with the 
Lee- Yang theorems and phase transitions are investigated for various choices of (3k- More 
complex iterative models of the partition function are proposed whose distribution of zeros 
are fractal sets. 

This paper also presents an alternative point of view for modeling fragmentation pro- 
cesses which is based on Markov process models. Markov process models give a picture 
of the underlying physical processes that lead to the cluster formation and break up. The 
relationship of this approach to that based on the canonical ensemble is discussed. 

Finally, some experimental data is investigated. Various choices for the quantity (3k are 
considered in our analysis. The statistics of the data are not sufficient to distinguish the 
various possible /Vs considered. 
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APPENDIX: GENERATING FUNCTIONS AND RECURRENCE RELATIONS 

In this section we consider a general procedure for generating the canonical partition 
function Qa{x,(J) from a generating function given by 

{OO 1 OO 

§£}||n"f = E^(§,})-tj ( A1 ) 
|| =00 ) A=i 

Here the g k s are arbitrary functions of k. Letting y = x J29kU k , expanding expy = 1 + y + 
y 2 /2 + . . . and collecting all terms with equal powers of u gives 

Qa(xJ) = £ M 3 {n)x m f[ (g k k\) nk (A2) 

U A (n) k=l 

where M 3 (n) = A\j Y\nk\k\ nk and consequently 

Qa{x, g) = E n .,;.-iw (A3) 

Thus Ph = g^ 1 . For example, 

Qi(x,g) = xg x 
Qi(x,g) = 2\g 2 x + gfx 2 

Q 3 {x, g) = 3\g 3 x + 3 9l (2\g 2 )x 2 + gfx 3 (A4) 
Q 4 (x,g) = (4\g 4 )x+\3(2\g 2 ) 2 + 4 gi (3\g 3 ) 



x 2 



+Qg 2 1 (2\g 2 )x 3 + gfx 4 

Using a procedure in Riordan |J, the following recurrence relationship is obtained 

Q A+1 (x,g) = Ixgi + V(s + l)g s+1 ^-\ Q A (x,g) (A5) 
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FIGURES 

FIG. 1. Permutations among particles. Left hand side of graph shows a group of particles and a 
permutation operator as it would act on the particles. Right hand side gives cluster interpretation 
of the same permutation. 

FIG. 2. The behavior of (m^) = k(nk) for the polymer model with A=50. (1) is the exact 
model given by fa = k\(k + 2)!/(2fe)!, (2) is the approximate model given by the Stirling limit 
fa = k 1 / 2 (k + l)(k + 2). 

FIG. 3. Zeros of the partition function Za(x) with q = — 1 

FIG. 4. Specific heat of a finite Bose gas for d=3. ksT\ = 8 MeV in this figure. 

FIG. 5. Zeros of the partition function Za(x) for the choice fa = k 2 scaled by \/A (i.e. root x 
is plotted at x/A). The cases A = 25, 50, 75, 100 are shown. 

FIG. 6. Zeros of the partition function Za(x) for the choice fa = k 3 scaled by 1/A (i.e. root x 
is plotted at x/A). The cases A = 25, 50, 75, 100 are shown. 

FIG. 7. The behavior of (m^) for the choice fa = k, A = 50 at various x. 

FIG. 8. The behavior of (m^) for the choice fa = 1, A = 50 at various x. 

FIG. 9. The behavior of (m^) for the choice fa = /c 1 / 2 , A = 50 at various x. 

FIG. 10. The behavior of (m^) for the choice fa = k 2 , A = 50 at various x. 

FIG. 11. The behavior of (m^) for the choice fa = k\, A = 50 at various x, with a comparison 
to the approximate model considered in the text. 
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FIG. 12. log(rifc) vs. k for best fits. Line (1) is for [5k = k. Lines (2), (3), (4), are for 
l/f3 k = a /k + (1 - a)/k T , with (2) T = 2, (3) r = 2.5, (4) r = 3. 
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TABLES 



TABLE I. Recursion relations for various fa 



TABLE II. Values of z ( A > in the large A limit 

Jl) _(2) _(3) (4) (5) 

r z A z A z A z A z A 

3/2 1.000 2.612 3.412 2.971 1.941 

2 1.000 1.645 1.353 0.742 0.305 
5/2 1.000 1.342 0.899 0.402 0.135 

3 1.000 1.202 0.723 0.289 0.0870 



TABLE HI. Markov Process Models 

Ajfe Hjk x fa 

a p (3/a 1 

a{jk) T faj + k) T j3 1 'a k T 

a !3{ j + k ) T [3 /a k\ T 
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TABLE IV. Restricted Markov process models 



X k fi k x /3k 

a (3 (3/a 1 

ak (3k f3 /a k 

a (3k (3/a k\ 

ak a fik T fi/a 



TABLE V. Fits of Experimental Data to various Models 

(3k x a t a 2 



n/a 78.58 

1.817 63.10 

1.5 65.61 

2 62.34 
2.5 61.55 

3 62.20 



k 0.296 n/a 

k T 1.895 n/a 

0.93745 

1.96263 0.0587 

2.61313 0.0747 

3.0803 0.0732 



a+(l- 






k 


a+(l- 


-a)k>-- T 




k 


o+(l- 


-a)k^ T 




k 

„M„1-T 
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